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Abstract 

We present Dirac-Brueckner-Hartree-Fock calculations for isospin asymmetric nuclear matter 
which are based on improved approximations schemes. The potential matrix elements have been 
adapted for isospin asymmetric nuclear matter in order to account for the proton-neutron mass 
splitting in a more consistent way. The proton properties are particularly sensitive to this adaption 
and its consequences, whereas the neutron properties remains almost unaffected in neutron rich 
matter. Although at present full Brueckner calculations are still too complex to apply to finite nu- 
clei, these relativistic Brueckner results can be used as a guidance to construct a density dependent 
relativistic mean field theory, which can be applied to finite nuclei. It is found that an accurate 
reproduction of the Dirac-Brueckner-Hartree-Fock equation of state requires a renormalization of 
these coupling functions. 

PACS numbers: 21.65.+f,21.60.-n,21.30.-x,21.30.Fe 

Keywords: Nuclear equation of state, isospin dependence, relativistic Brueckner approach, density dependent 
relativistic mean field theory 
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I. INTRODUCTION 



The investigation of asymmetric matter is of importance for astrophysical and nuclear 
structure studies. In the field of astrophysics this investigation is important for the physics of 

il composition and cooling 



6J, and some other topics. 



supernova explosions [lj and of neutron stars [2j, e.g. the chemica 
mechanism of protoneutron stars 0, 4], mass-radius correlations 5, 
In the field of nuclear structure the investigation of isospin asymmetric matter is of interest 
in study of neutron-rich nuclei |7|. This isovector dependence of the nuclear force can be 
investigated in the heavy ion experiments 8] . However, the data for neutron-rich nuclei were 
rather scarce in the past. This situation is changing with the forthcoming new generation of 
radioactive beam facilities, e.g. the future GSI facility FAIR in Germany, the Rare Isotope 
Accelerator planned in the United States of America or SPIRAL2 at GANIL/France, which 
will produce large amounts of new data. 

Models which make predictions on the nuclear equation of state (EoS) can roughly be 
divided into three classes: Phenomenological density functionals, effective field theory (EFT) 
approaches, and ab initio approaches. Phenomenological density functionals are based on 
effective density dependent interactions such as Gogny or Skyrme forces or relativistic 
mean field (RMF) models 



10l | with usually more than six and less than 15 parameters. The 



effective field theory approaches lead to a more systematic expansion of the EoS in powers of 
density, respectively the Fermi momentum hp. The advantage of EFT is the small number 
of free parameters and a correspondingly higher predictive power. Ab initio approaches are 
based on high precision free space nucleon-nucleon interactions and the nuclear many-body 



problem is treated microscopically. Predictions 
free. Examples are variational calculations 



or the nuclear EoS are essentially parameter 



111 ]. Brueckner- 



14 



15 



artree-Fock 



16 



la nana 



(BHF) 



or 



20J calculations 



2l|. 



relativistic Dirac-Brueckner-Hartree-Fock (DBHJ 
and Greens functions Monte-Carlo approaches 

Many-body calculations, on the other hand, have to rely on the summation of relevant di- 
agram classes and are still too involved for systematic applications to finite nuclei. However, 
these results can be used as a guidance for the construction of a "semi" -phenomenological 



density functional. Examples are e.g. Gogny forces [22 



dependent relativistic mean field (DDRMF) theory 
results. 



23, 



derived from G-matrices or density 



241 ]. which can be based on DBHF 
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The theoretical predictions for the isospin dependence of nuclear interactions are still 
diverse. RMF theory can not describe the complex nonlinear behavior of the DBHF and 
BHF binding energy at densities near p = 0. Furthermore, the symmetry energy in rela- 
tivistic DBHF calculations is found to be significantly stiffer than in non-relativistic BHF 
approaches [25j , in particular at high densities. The BHF calculations 26( predict a proton- 
neutron mass splitting of ni* NRn > m* NRp in neutron dominated nuclear matter. In con- 
trast, RMF theory with the scalar isovector 5-meson included predict the opposite behavior, 



m D,n < m u, P 



*D,P 



271 ] . The various Skyrme forces give opposite predictions for the neutron- 
proton mass splitting. Relativistic ab initio calculations based on realistic nucleon-nucleon 
interactions, as for instance the DBHF approach, are the proper tool to answer these ques- 
tions. 

In this work we describe asymmetric nuclear matter in the framework of the relativistic 
DBHF approach based o n p rojection techniques using the Bonn potential and their bare 
NN matrix elements V 2a] . Furthermore, the optimal representation scheme for the T- 
matrix, the subtracted T matrix representation, is applied. This scheme has previously 

18, 



been applied to asymmetric nuclear matter in refs v 



work we go beyond the approach used in 



18 



[l9 . 2oj|. However, in the present 



in the sense that we improve at a 



couple of approximations. To be more precise, the Bonn potential has now been adapted 
for asymmetric nuclear matter. 

In the solution of Bethe-Salpeter (BS) equation we abandon the approximation of an 
averaged np mass in the np channel and distinguish explicitly between the different isospin 
dependent matrix elements. As a consequences, the potential and T-matrix are evaluated 
in terms of six independent helicity or covariant amplitudes instead of five [151 ]. which are 
sufficient in the case of an averaged np mass. 

The plan of this paper is as follows. The relativistic DBHF approach with emphasis on 
the treatment of the nn, pp, and np channels is treated in sect. [Ill Results are presented in 
sect. IIII1 Furthermore, the relation between DBHF results and the RMF theory is discussed 
in sect. IIVBI Finally, we end with a summary and a conclusion in sect. [V] 
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II. DBHF APPROACH IN ISOSPIN ASYMMETRIC NUCLEAR MATTER 



In this section the relativistic Brueckner approach is discussed. First a general overview 
is given, followed by a more detailed discussion of the modifications, which are necessary to 
account properly for the proton-neutron mass splitting and the isospin dependence of the 
corresponding matrix elements. 

In the relativistic DBHF approach a nucleon inside nuclear matter is regarded as a dressed 
particle as a consequence of its interaction with the surrounding nucleons. This interaction 
of the nucleons is treated in the ladder approximation of the relativistic BS equation 

T = V + i J VQGGT, (1) 

where T denotes the T-matrix, V the bare nucleon-nucleon, Q the Pauli operator, and G 
the Green's function of an intermediate off-shell nucleon. This Green's function G which 
describes the propagation of dressed nucleons in the medium fulfills the Dyson equation 

G = G + G SG, (2) 

where Go denotes the free nucleon propagator and S the self-energy. In the Hartree-Fock 
approximation this self-energy is given by 

Z = ~iJ [Tr[GT] - GT). (3) 

F 

The eqs. (EE])-© are strongly coupled. Therefore, this set of equations represents a self- 
consistency problem and has to be iterated until convergence is reached. 

The structure of the self-energy follows from the requirement of translational and rota- 
tional invariance, hermiticity, parity conservation, and time reversal invariance. The most 
general form of the Lorentz structure of the self-energy in the nuclear matter rest frame is 
given by 

E(k, k F ) = E B (k, k F ) - 7 E Q (k, k F ) + 7 • k E v (k, k F ), (4) 

where the S s , S Q , and S v components are Lorentz scalar functions which depend on the 
Lorentz invariants k 2 ,k ■ j and j 2 , with the baryon current. Therefore, these Lorentz in- 
variants can be expressed in terms of fco, |k| and k F , where &f denotes the Fermi momentum. 
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The different components of the self-energy are determined by taking the respective traces 



29. 



30| 



£ B = ^r[£], £ =-Ur[ 7 o£], S v = — ^tr [7 • kS] . (5) 



The presence of the medium influences the masses and momenta of the nucleons inside 
nuclear matter. These effective masses and effective momenta of the nucleons can be written 

as 

m*(k,k F ) = M + $teZ s (k,k F ), k*^ = hP + 3te£"(ife, k F ). (6) 
By the introduction of reduced quantities, one has the reduced effective mass 

m*(k, k F ) = m*(k, fcp)/ (1 + E y (k, k F )) . (7) 
and the reduced kinetic momentum 

k; = k;/(i + z v (k,k F )) , (8) 

Hence, The Dirac equation written in terms of these reduced effective masses and momenta 
has the form 

py^fe*" - m* (k, k F )}u{k, k F ) = 0. (9) 

To simplify the self-consistency scheme we will work in the quasi-particle approximation, 

re self-energy SmS will be neglected. In addition, the "reference 



i.e. the imaginary part of t 
spectrum approximation" |3l|] is applied in the iteration procedure, i.e. the effective mass 
of the nucleon is assumed to be entirely density dependent (|k| = kp). However, in general 
the reduced effective mass is density and momentum dependent. Therefore, this method 
implies that the self-energy itself is only weakly momentum dependent. At the end of the 
calculation one has of course to verify the consistency of the assumption ~ S(|k| = k F ) 
with the result of the iteration procedure. 

The solution of the Dirac equation in eq. provides the positive-energy in-medium 
nucleon spinor 



'E*(k)+m* F / 1 



^ M = \l v 4* F 2A | k , I Xa, (io) 

* £*(k)+mj 



where E*(k) = yk 2 + m*p denotes the reduced effective energy and xx a two-component 
Pauli spinor with A = ±| 44]. The normalization of the Dirac spinor is thereby chosen as 
u x (k, k F )u\(k, k F ) = 1. 



A. nn and pp channel 

It is convenient to reduce a four- dimensional BS integral equation, eq. |T|, to a three- 
dimensional one to solve the scattering problem of two nucleons in the nuclear medium. 
Therefore, the two-particle propagator iGG in the BS equation has to be replaced by the 
effective Thompson propagator. The Thompson propagator implies that the time-like com- 
ponent of the momentum transfer in V and T is set equal to zero. Hence, the Thompson 
propagator restricts the exchanged energy transfer by 5(k°) to zero. In addition, the Thomp- 
son propagator projects the intermediate nucleons onto positive energy states. Thus, in the 
two-particle center of mass (cm.) frame, which is the natural frame for studying two-particle 
scattering processes, the Thompson equation can be written as 



T(p,q,x)U. = V{p,q) (11) 

f d * k a TV ia m F Q(k,x) 

+ J (2vr)3 VlP ' kj £*2(k)2£*(q) -2E*(k) + le ( 

where q = (q x — q 2 )/2 is the relative three-momentum of the initial state and k and p are 
the relative three-momenta of the intermediate and the final states, respectively. 

The Thompson equation (fTTT) can be solved applying standard techniques, which are 
outlined in detail by Erkelenz 32| . To determine the self-energy only positive-energy states 
are needed. Therefore, it is more convenient to apply the Dirac nucleon propagator [29! ]. 

G D (k, k F ) = [7^*" + m*{k, k F )]27Ti5{k* 2 - m* 2 {k, k F ))Q{k*°)Q{k F - |k|)., (12) 

instead of the full nucleon propagator. Due to the 9-functions in eq. ( TT2|) only positive 
energy nucleons are allowed in the intermediate scattering states. In this way, one avoids 
the delicate problem of infinities in the theory which generally will occur if one includes 



contributions from negative energy nucleons in the Dirac sea 



15 



m. 



In the on-shell case for identical particles only five of the sixteen helicity matrix elements 



are independent which follows from general symmetries 32j. After a partial wave projection 



onto the | J MLS >-states the Thomas equation reduces to a set of one-dimensional integral 



equations over the relative momentum |k|. Furthermore, it decouples into three subsystems 
of integral equations: the uncoupled spin singlet, the uncoupled spin triplet, and the cou- 

integral 



29|. Due 



pled triplet states (appendix [B]). To achieve this reduction to the one-dimensiona . 
equations the Pauli operator Q is replaced by an angle-averaged Pauli operator Q 
to deformation of the Fermi sphere to a Fermi ellipsoid in the two-nucleon cm. frame, Q is 
evaluated for such a Fermi ellipsoid: 

Ik I <r h ■ 

Q = { " E 1L E ' F for k mm < |k| < k max (13) 

h <r Ikl 




with k min = y/kp — u 2 Ep, k max = ^{uEp + k F ), and u = |u|. The partially decoupled set 
of the one-dimensional integral equations are solved by the matrix inversions techniques of 
Haftel and Tabakin [3^]. 

Due the anti-symmetry of these two-nucleon states the total isospin of the two-nucleon 
system (I = 0, 1) can be restored by the standard selection rule 

(_l) L +s+i = _i. (14) 

The five independent partial wave amplitudes in the helicity representation are obtained 
from the five independent on-shell amplitudes in the \JMLS >-representation [32|]. After 
the summation over the total angular momentum one has the five on-shell plane-wave helicity 
matrix elements 

< P A' 1 A' 2 |T I (x)|qA 1 A 2 >= J2 (^r) 4a(*) < A' 1 A;|T J ' I (p,q,x)|A 1 A 2 >, 

(15) 

where 6 is the scattering angle between q and p with |p| = |q|. Furthermore, one has 
A = Aj — A 2 and A' = A x — A 2 . The reduced rotation matrices d 3 , (9) are those defined by 
Rose y. 



Since we determine the T-matrix elements in the two-particle cm. frame, a representa- 
tion with covariant operators and Lorentz invariant amplitudes in Dirac space is the most 
convenient way to Lorentz-transform the T-matrix from the two-particle cm. frame into the 
nuclear matter rest frame 29|. Some freedom in the choice of this representation exists, be- 



cause pseudoscalar (ps) and pseudovector (pv) components can not uniquely be disentangled 



for on-shell scattering. This ambiguity is minimized by separating the leading order, i.e. the 
single-meson exchange, from the full T-matrix. Therefore, the contributions stemming from 
the single-7r and-n exchange are given in the complete pv representation. For the remaining 
part of the T-matrix, the ps representation is chosen. 

Taking the single nucleon momentum k = (0, 0, k) along the z-axis, then we have for the 
nn and pp channel contributions for the self-energy components in the ps representation 
scheme 

1 r kp i <7 3 n m* 



S?(k) = \ fj 3 ^["4if + Ff-2F»-2F»- if], (17) 



and 



1 fkFj ^/3 n n u 



1 ftei d 3 q q-k 
q,j 

where % = j = n or % = j = p, respectively. In the complete pv representation the nn and 
pp channel contributions to the self-energy components are given by 

,*2 1 ^*2 



S/(k) = 4 X (^F^ W - * + 4 ^ + (rnf + ro;)a ^ 



(19) 



Lo(k) " + 4i (2vr)3^ 2 <? A i^K + m*) 2 ^(20) 



and 



1 d 3 q q-k 



2k* (mf - k*^q*^) - q z (m* 2 - m* 2 
q z (m* + m J 



1 1 ij 



,*\2 ^pv 



9U ( 21 ) 



where i = j = n or i — j — p, respectively. 

In short, the complete pv representation is applied for V^ j7J and the representation is 
used for the T Sub = T — V4, T; to get the most favorable representation scheme, the subtracted 
T-matrix representation scheme. 
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B. np channel 



Since in isospin asymmetric nuclear matter one has to deal with two distinct nucleons 
states in the np channel, this channel is more complicated than the nn and pp channel. 
Working with two distinct nucleons has consequences for the Thompson equation, the Pauli 
blocking operator, and the number of independent helicity matrix elements. 

First, the Bonn potential 



281 ] has to be made suitable to treat distinct particles in the 
medium. An important difference is that the neutrons and protons have unequal effective 
masses. These distinct effective masses have to be accounted for, in particular in the eval- 
uation of the potential matrix elements. The resulting one boson exchange (OBE) matrix 
elements can be found in appendix [A] 

Second, the two-particle propagator iGiGj in the BS equation has to be replaced by 
the Thompson propagator for the np channel. The effective Thompson propagator for this 
channel is given by 

g np - l^n^p - — — — — — — , \2.2.) 

K E p Vs* - E* - E* + le 



or iden- 



where \fs* is the invariant mass. 

In contrast to the five independent helicity matrix elements in the on-shell case 



tical particles, in the np channel six helicity matrix elements are independent [35j|. After 
the partial wave projection onto the \ JMLS >-states, using an average direct-exchange 
contribution in the potential this time the Thompson equation partially decouples into two 
subsystems of one-dimensional integral equations: the coupled spin singlet-triplet states and 
the coupled triplet states (appendix [B]). To achieve the reduction to the one-dimensional 
integral equations the Pauli operator Q has to be replaced by an angle-averaged Pauli oper- 



ator Q 



29j . However, the Pauli operator Q for the np channel has to be modified compared 



to the one in the nn and pp channel, since it has to be evaluated for Fermi ellipsoids with 
different sizes. The result for the angle-averaged Pauli operator for the np channel Q np with 
a neutron excess is 



Q 



up 



Q(-yuE Fn - 7/cFn) |k| < k min 

l/2[cos(0 p ) - cos(6> n )]9(6» n - 9 p ) for k min < |k| < k max (23) 

1 krn a nr. *^ k 
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with k min = "i\uE Fn — k Fn \, k m ax = l(uE Fn + k Fn ) 



and 



, arccosf^f^) for |^^| < 1 
6 P = { V 7l k IM J J 1 Tl k H u l 1 " , (24) 

otherwise 



e n = { v 7M]u| ; j i 7 |k|| U | i - _ (25) 

7r otherwise 



Due to the additional independent helicity matrix element, we will have a sixth indepen- 



dent covariant in the T-matrix representation 35(. However, the problem is that we need 



to have a decomposition that reduces to the one used in the symmetric case. The general 



Lorentz representation of the nine invariants given in [35| fulfill this requirement. Leaving 
out the three redundant invariants in our case, the additionally constructed covariant is 
defined as 



Tl> dir (\p\,e,x) = li^cipi, *)((tm)2 • Of - (tm)i - Q 



1 F 6 I (|p|^,x)( 72 -k + 7i-k), (26) 



2 

with = (pi + q' i )' 1 /2m for % = 1, 2 and k = p x + q' x = — (p 2 + q' 2 ) in the cm frame. The 
same sixth covariant is used in ref . [ljj] , while the exchange sixth amplitude given in ref . 3] 
does not contribute in the self energy components. Therefore, one gets an additional term 
in the np channel contribution to the neutron self energy components 

s n P(k) _ 1 [*» d»q m; k* % -m? 



1 f kFp d 3 a Et - E* 



and 

compared to the nn and pp channel. For the proton a similar additional term arises, where 
neutrons and protons are interchanged in eqs. (p7|) - (!29|) . In symmetric nuclear matter with 
equal effective masses for neutrons and protons, the coefficient of this sixth independent 

10 



amplitude vanishes, i.e. the familiar representation scheme with the five linearly independent 
covariants is obtained, as expected. 

Finally, the total neutron and proton self energies including all channels can be written 

as 



j-m ^itm _|_ Y^ n P ■ VJPP _|_ YjP n 



(30) 



respectively. 



III. RESULTS AND DISCUSSION 
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FIG. 1: Binding energy as a function of the baryon density. 



In fig.[T]we present the results for the equation of state for various values of the asymmetry 
parameter (3 = (n n —n p )/riB in the framework of the DBHF approach with a sixth indepen- 
dent amplitude in the np channel using the Bonn A potential. The applied representation is 
the optimal representation so far, the subtracted T-matrix representation. The two extreme 
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he symmet 



16 



18] ■ The 



cases are symmetric nuclear matter (/3 = 0.0) and neutron matter (f3 = 1.0). 
ric nuclear matter results and neutron matter results agree with those of refs. 
binding energy curves for intermediate values of f3 lie between these two extreme curves and 
are slightly higher than in ref. [18] . In addition to that, the binding energy 



E{n B , (3) = E{n B ) + E sym {n B )f3 2 + 0(p 4 ) 

shows a nearly quadratic dependence on the asymmetry parameter (3 as expected. 

150 i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 



(31) 
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FIG. 2: Comparison o:' 
a nonrelativistic BHF 



several EoSs from ab initio calculations, i.e. the 
I2I ] (dashed-dotted) and a variational calculation 



aresent approach (solid), 
11[ (dotted). 



Fig. [2] compares our prediction for the binding energy to the ones of other microscopic 



many-body approaches, the variational calculations from |TJJ and the nonrelativistic BHF 
calculation from [l2|. at symmetric nuclear matter (below zero) and pure neutron matter 
(above zero). The variational calculation is based on the high precision phenomenological 



Argonne Vi 8 [36j two-nucleon interaction and includes UIX three-body forces [37| as well as 



12 



relativistic boost correction denoted by Sv 11] . Also the nonrelativistic BHF calculation 



12 



is based on the phenomenological Argonne Vis |36j. Furthermore, it includes a microscopic 
three-body force deduced from the meson-exchange current approach [l^ . 

The first observation that becomes evident from fig. [2] is that in nuclear matter both, 
the BHF and DBHF, calculations lead to more binding than the variational calculation. 
However, in all three cases the EoS of nuclear matter can be characterized as "soft', at 
least at moderate densities up to about three times saturation density. The prediction of a 
soft EoS is the general outcome of a microscopic many-body calculation. Recent Quantum 
Monte Carlo calculations for symmetric nuclear matter [38( show the same tendency. It 
should be noticed that this observation is supported by corresponding observables extracted 
from heavy ion reactions, where supranormal densities up to about three times saturation 
density are probed. Heavy ion data for tranverse flow [39( or from kaon production 40] 
support the picture of a soft EoS in symmetric nuclear matter. 

In neutron matter the variational calculations are less stiff, in particular at high-density 
neutron matter, than our DBHF calculations, whereas the nonrelativistic BHF calculation 
lies in between these two approaches. However, up to 1.5 times saturation density for neutron 
matter and symmetric nuclear matter the three approaches show a quite reasonable agree- 
ment. This fact indicates that these are the density ranges which are at present reasonable 
well controlled by state-of-the-art many-body calculations. 

The high density behavior of the EoS, in particular that of the neutron matter EoS, can 
be constructed by astrophysical observables Q|. The recent observation of the at present 
heaviest compact star, a binary pulsar of 2.1 ± O.2M (la level) 4jJ rules out very soft 
neutron matter EoSs. However, all three EoSs shown in fig. [2] fulfill this constraint since 
they yield maximum neutron star masses between 2.2 2.3M . 

In fig. [3] the neutron and proton optical potentials are plotted as a function of the mo- 
mentum k = |k| for various values of the asymmetry parameter (3 = (n n — n p )/riB at a 
fixed nuclear density of % = 0.166 fm~ 3 . The depth of neutron optical potential decreases 
with increasing asymmetry, whereas the depth of the proton optical potential shows the 
opposite behavior. Furthermore, the steepness of the neutron optical potential decreases 
with increasing asymmetry parameter (3, whereas the opposite behavior is found in the pro- 
ton case. Compared to ref. {20] the neutron optical potential remains almost unaltered. In 
contrast, the proton optical potential lies a bit lower and is somewhat steeper as compared 
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FIG. 3: The neutron and proton optical potential in neutron rich matter as a function of the 
momentum k = Ikl. 



to ref. 

The isovector optical potential Ui SO = ^-fi strongly depends on density and momentum. 
This optical potential in neutron-rich matter initially stays constant and then decreases 
strongly with increasing momentum. Furthermore, the isovector optical potential is almost 
independent of the asymmetry parameter (3. This behavior can also be observed in refs. 18, 



19 



20| | . Since the proton optical potential lies a bit lower, the isovector optical potential at 



k = is slightly higher than in refs. is[ [l9j, |20j. However, the optical isovector potential 
at nuclear density ub = 0.166 fm -3 at k = is still in good agreement with the empirical 
value of 22 - 34 MeV [42 1. 



An interesting issue is the proton- neutron mass splitting in neutron-rich matter, which has 
in detail been discussed in refs. 



19 



20j. One should keep in mind that different definitions 



of the effective mass exist, which are often compared and sometimes even mixed up: the 
nonrelativistic mass and the relativistic Dirac mass. In fig. [4] the nonrelativistic and Dirac 
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FIG. 4: Neutron and proton effective mass as a function of the momentum k = |k| in neutron 
matter at fixed nuclear density % = 0.166 fm -3 . In addition, the effective mass in symmetric 
nuclear matter is given. 



effective mass of the neutron and proton are compared for (3 = 1, i.e. neutron matter. Our 
DBHF calculations based on projection techniques predict a mass splitting of m* Dn < tn* D 
in neutron-rich matter. However, the nonrelativistic mass derived from our DBHF approach 
shows the opposite behavior. This opposite behavior to the relativistic Dirac mass, i.e. 
rn% Rn > m* NRp , is in agreement with the results from nonrelativistic BHF calculations 
261 ] . This difference between the Dirac mass splitting and the nonrelativistic mass splitting 
is not surprising, since these masses are based on completely different physical concepts. 
The relativistic Dirac mass is defined through the scalar part of the nucleon self-energy in 
the Dirac field equation which is absorbed into the effective mass ([6]). On the other hand, 
the nonrelativistic mass parameterizes the momentum dependence of the single particle 
potential. 

In this context we want to note that, in contrast to the non-relativisitc mass m* NR , the 
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FIG. 5: Neutron and proton effective mass as a function of the density in neutron matter. In 
addition, the effective mass in symmetric nuclear matter is given. These results are compared to 



DBHF calculations with average mass assumption in the np channel (5 cov.) of refs. 
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20]. 



momentum dependence of the Dirac mass m* D is smooth and still moderate. This fact is 
important to justify the reference spectrum approximation, i.e. the usage of an momentum 
independent effecive Dirac mass fh* F for the evaluation of the in- medium spinor basis fllOl) . 
the Thompson propagator and the potential matrix elements (see Appendix). 

In fig. [5] the neutron and proton effective Dirac masses are plotted as a function of the 
baryon density % for pure neutron matter. Of course, a strong density dependence can be 
observed. In addition, one has a Dirac mass splitting of m* D n < m* D in the whole density 
range. RMF field theories with the isovector p and S mesons included predict the same Dirac 
mass splitting. When only the p-meson is included, the RMF theory predicts equal masses. 
Hence, the S meson is responsible for the mass splitting in RMF theory. 

Furthermore, in fig. [5] our results for the neutron a nd p roton effective Dirac mass in pure 
neutron matter are compared to those from refs. H, 2o| , where only 5 covariants were 
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used in the np channel. The neutron effective mass remains practically unaffected, whereas 
the proton mass experiences a sizable reduction. These results are easy to understand. The 
neutron self energy consists of a nn and a np part. Hence, the nn part becomes dominant 
for a vanishing proton fraction. The proton self energy consists of a pp and a np part. In the 
limit of a vanishing proton fraction, the np interaction becomes dominant. Therefore, the 
proton properties, e.g. the proton effective mass, are especially sensitive for the treatment 
of the np channel in neutron rich matter. 



IV. RELATION TO RELATIVISTIC MEAN FIELD THEORY. 



A. DBHF self energy components 

At present full Brueckner calculation are still too complex to allow an application to fi- 
nite nuclei. However, within the framework of density dependent mean field theory effective 
density depend coupling functions can be obtained from the Brueckner self-energy compo- 
nents. Such coupling functions parameterize the correlations of the T-matrix in a handable 



way and can be applied to finite nuclei within the framework of DDRMF theory 23]. In 
contrast to standard RMF models, the meson-baryon vertices are density dependent. As a 
consequence, rearrangement contributions in the baryon field equations occur. These rear- 
rangement contributions should be taken into account and are essential to satisfy energy- 
momentum conservation and thermodynamic consistency in this density dependent mean 
field theory. 

In order to properly parameterize the isospin dependence of the self-energy components, 
the coupling functions must be based on four different channels: scalar isoscalar, vector 
isoscalar, scalar isovector, and vector isovector channel. In RMF theory these channels 
correspond to phenomenological exchange bosons, i.e. the er, uo, 5, and p mesons. The 
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effective coupling constants are then given by 



m a ) 2 n s 
g<j(nB,(3)\ 2 _ 1 S 0)P (A; Fp ) + S 0in (A; F „) 



J 2 n B 
9s(nB, P) \ 2 _ 1 S SiP (A; F p) - E B , n (k Fn ) 



m$ J 2 n s3 
( 9 P (n B ,P)\ 2 _ 1 S 0jP (A; F p) - S 0in (/c Fn ) 



\ m p J 2 n 3 

with n s = ra sp + n sn , n B = n p + n n , n s3 = n sp - n sn , and n 3 = n p - n n , where 



(32) 
(33) 
(34) 
(35) 



2 f kFt „ m? 



n si = - — - / dk , 1 = (36) 
Jo y/m? + k 2 1 ^ 



and 

-k Fi ^,3 



^ = 77T^ / dk =^H (37) 



(2tt) 3 7 3tt 2 

are, respectively, the scalar and vector density of particle i(= n,p). The results for the 
isoscalar and isovector coupling constants are plotted in figs. [6] and [7J respectively. The 
strength of the isoscalar coupling functions decreases as the density increases. At low den- 
sities, both the scalar g a and the vector isoscalar coupling g u show a strong decrease with 
increasing density. However, at higher densities the vector coupling stays more or less con- 
stant. The strength in the isovector channel is small compared to that in the isoscalar 

n 

channel. Furthermore, compared to ref. [18J the dependence on the proton fraction for 
isovector strength is strongly reduced. 

We can directly use these density dependent coupling functions in a RMF theory for 
infinite nuclear matter. For practical purposes and to keep the DDRMF functional as 
simple as possible we ignore in the following the weak isospin dependence and assume only 
a density dependence in the coupling functions of eqs. (132]) -( 135]) . In fig. [8]the binding energy 
determined from this RMF theory is compared to our DBHF results for neutron matter and 
symmetric nuclear matter. The binding energy in RMF theory is given by 

2 



ri'i 



i=n,p 

\ 2 / r, ( m ^\ \ 2 

2 



+ r-M) 2 ^ + (^) 2 n 3 



(38) 
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FIG. 6: The isoscalar scalar (g a ) and vector (g^) effective coupling functions are plotted as a 
function of the baryon density for different values of the asymmetry parameter j3. 



where the density dependent couplings 



9a(n B ) 



m r 



,a e {a,u,5,p} 



(39) 



are obtained from eqs. (I32j) - (l35l) using the data for f3 = 0.2. No rearrangement terms are 



present in eq. 



since rearrangement contributions do not contribute at the level of the 



binding energy 23]. In RMF theory the integral for the kinetic energy can be evaluated and 



leads to the analytical expression 
2 



(2k) 



E / d3kE ^ = E 

=n,p Fi i=n,p 



3 p 1 * 

-E Fi ni + -17^71, 



(40) 



with the Fermi energy E Fi = ■sj k 2 Fi + m* 2 . The effective mass contains the contributions of 
the two scalar mesons. Through the different coupling to the isovector 5-meson it accounts 
for the proton-neutron mass splitting, i.e. 



m 



n/p 



M — 



m s 



n S 3 



(41) 
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FIG. 7: The isovector scalar (g$) and vector (g p ) effective coupling functions are plottted as a 
function of the baryon density for different values of the asymmetry parameter (3. 



Comparing the original DBHF EoS in Fig (JSj) for the DDRMF EoS based on the parameter- 
ization (|32|) - (|35|) . one observes clear deviations of the two approaches, both for symmetric 
as well as for neutron matter. This suggests that the density dependent coupling functions 
should be extracted more carefully as has been done in the "naive" definition (I32l) -( l35l) . 
With other words, an accurate reproduction of the DBHF EoS requires a renormalization of 
the coupling functions which includes the contributions from Fock terms in a more consistent 
way. 



B. Renormalized self energy components 



The fact that renormalization is required when DBHF results are mapped on RMF theory 
can easily be seen from the DBHF binding energy, 

2 x - 



E/A 



(27T 



|3 A^/ 

i=n,p 



(42) 
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FIG. 8: The DBHF EoS is compared to DDRMF EoS and the renormalized DDRMF EoS. 



The two essential differences between DBHF and RMF concerning the structure of the self- 
energy, respectively the mean field, are firstly, that the DBHF self-energies carry an explicit 
momentum dependence and, secondly, the appearance of a nonvanishing spatial contribution 
£y, see eqs. (JI])-©. Both features should be taken into account as accurate as possible when 
DBHF results are parameterized in terms of RMF theory. The £y component originates 
from Fock exchange contributions which are not present in RMF theory. For an accurate 
reproduction of the DBHF energy functional the spatial £y component has to be included 
in proper way. Firstly, the £y component can be absorbed into the effective mass according 
to eq. ([7]) and this reduced effective mass has to be identified with RMF effective mass, i.e. 

= M + S„, t (k Fl ) =M y^DDRM F (43) 

This leads to the renormalized scalar self energy component 

y DDRMF _ E s ,t(fc_Ft) ~ MH^^kpi) . ... 

l + Ev,i(A*i) ' 1 1 
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However, the DBHF energy functional of eq. (142]) has additional terms compared to the 
DDRMF energy functional of eq. fl38l) . In the same way, however, then using the energy 
density instead of the effective mass, the following expression for the normalized vector self 
energy component is obtained 

^v,i(k Fi )[3E Fi ni + fh*n Sji ) 



^DDRMF 
J 0,i 



Fi 



4n,- 



(45) 



These renormalized self energy components are now inserted into eqs. (1321) - (1551) to obtain 
the renormalized density dependent coupling functions. By this procedure all terms which 
contribute to the DBHF energy functional are taken into account in the correspondingly 
constructed DDRMF functional. However, the explicit momentum dependence in eq. ( 142]) 
can not so easily be transferred to the RMF theory which leads still to slight deviation of the 
corresponding energy functionals. A possibility would be to perform a Taylor expansion of 
the self-energy components in terms of the momentum 24j. Since the intrinsic momentum 



dependence of the DBHF self-energy components is generally weak we neglect such 

additional correction terms. The new renormalized isoscalar density dependent coupling 
functions are reduced by an amount of 15-20 MeV compared to the corresponding nonrenor- 
malized coupling functions in fig. [6], but the qualitative behavior is very similar. The small 
isospin dependence of the isoscalar strength is almost insensitive to the renormalization. The 
renormalized isovector p and 5 meson coupling functions are shown in fig. The density 
dependence of the renormalized isovector coupling functions is similar to those of fig. 
However, the /3-dependence is now more pronounced. 

These renormalized density coupling functions are applied in RMF theory. Again the 
weak isospin dependence of these coupling functions is ignored. Therefore, we use again 
the data for j3 = 0.2 which corresponds approximately to the asymmetry in an Au nucleus. 
From table (TJ one can see that the saturation density is shifted to lower densities and the 



model 


n^ffm -3 ] 




E/A [MeV] 


DDRMF 


0.143 


1.28 


-14.22 


renormalized DDRMF 


0.168 


1.36 


-15.18 


DBHF 


0.181 


1.39 


-16.15 



TABLE I: Saturation properties of the DBHF model and the correspondingly constructed rela- 
tivistic mean field functionals. 
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FIG. 9: The renormalized isovector scalar (g$) and vector (g p ) effective coupling functions are 
plotted as a function of the baryon density for different values of the asymmetry parameter (3. 



binding energy of the saturation point is weaker in renormalized DDRMF theory compared 
to the original DBHF results. Without renormalization the deviations for the saturation 
density and the binding are even stronger. In fig. [8] the binding energy of the renormalized 
DDRMF theory is shown for neutron matter and symmetric nuclear matter. The results 
are in a fairly good agreement with the DBHF results and much better than the results 
without renormalization. Although the renormalized isovector coupling functions g 2 and g\ 
show a stronger /3-dependence the assumption of only density dependent couplings is still 
a good approximation. Extracting those coupling functions at the representative value of 
(3 = 0.2 both, the symmetric and the neutron matter EoS, are reproduced with fairly good 
accuracy. However, as can be seen from Table [J the mapping of DBHF onto the RMF 
functional leads generally to a shift of the saturation point towards lower densities and a 
slightly smaller binding energy. This feature which is mainly due to the neglection of the 
intrinsic momentum dependence of the DBHF self-energy has also been observed in previous 
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works when a similar procedure was applied 23|, |24j. Although the most realistic, Bonn A 
leads in DBHF to a slightly too large saturation density 14j, [16( and therefore such a shift 
of the saturation point is in favor of the DDRMF functional when applied to finite nuclei. 



V. CONCLUSION 

In summary, we present calculations of isospin asymmetric nuclear matter in a relativistic 
DBHF framework based on projection techniques. The approximation scheme for the treat- 
ment of isospin asymmetric nuclear matter has been improved. First of all, the application 
of the Bonn potential - factually the Bonn A potential has been used throughout this work - 
has been modified in order to distinguish between different proton and neutron mass by the 
evaluation of the potential matrix elements. The modification is essential in the np channel 
when in-medium matrix elements are evaluated. 

Secondly, the T-matrix can be represented by a set of six linearly independent Lorentz 
invariants in the np channel. This sixth covariant has been chosen as proposed in [jjj] . How- 
ever, in contrast to [15| we apply still the improved decomposition scheme of the T-matrix 



(subtracted T-matrix representation) 



16. 
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13, |2Cj which minimizes on-shell ambiguities 



in the determination of the self-energy components. 

It is found that the proton properties are, in particular, sensitive to the consequences 
of the adaption of the Bonn potential for isospin asymmetric nuclear matter and the intro- 
duction of a sixth covariant. The proton optical potential lies a bit lower and is steeper as 



compared to ref. 20j, whereas the neutron optical potential is almost unaltered. Further- 
more, the neutron effective mass remains practically unaffected, whereas the proton mass 
experiences a sizable reduction. The reason is that in neutron rich matter proton properties 
depend much stronger on contributions from the np channel than neutron properties. 

The main properties of isospin asymmetric nuclear matter remain, however, unchanged. 
The binding energy shows the expected quadratic dependence on the asymmetry parameter. 
Also the depth and the steepness of neutron optical potential decreases with increasing 
asymmetry, whereas the depth and steepness of the proton optical potential still shows the 
opposite behavior. A strong density and momentum dependence can again be observed 
for isovector optical potential. In addition, the isovector optical potential remains almost 
independent of the asymmetry parameter (3. Our DBHF calculations based on projection 
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techniques predict a mass splitting of m* D n < m* D p in neutron-rich matter as expected. The 
nonrelativistic mass derived from our DBHF approach still shows the opposite behavior, 
which is in agreement with the results from nonrelativistic BHF calculations 26]. 

At present full Brueckner calculations are still too involved for systematic applications in 
finite nuclei. However, the density dependent mean field effective coupling functions, which 
are obtained from the Brueckner self-energy components, parameterize the correlations of 
the T-matrix in a handable way. Therefore, these coupling functions can be applied to finite 
nuclei within the framework of DDRMF theory 23) . Doing so, a "naive" parameterization of 
the DBHF results in terms of a density dependent relativistic mean field functional leads to a 
poor reproduction of the original EoS. The reason are contributions from Fock-terms which 
are not present at the mean field level and which have to be incorporated in an effective 
way. This leads to a renormalization procedure of the coupling functions which absorbs the 
contributions from the Fock terms. However, the results with renormalization are in a quite 
good agreement with the DBHF results. With other words, an accurate reproduction of the 
DBHF EoS requires a renormalization of the coupling functions. 
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APPENDIX A: POTENTIAL MATRIX ELEMENTS. 



In this appendix we give the potential matrix elements for scalar, pseudovector, and vector 
mesons. Overall factors in front are omitted. In contrast to the normal expressions which are 
used in the Bonn codes [43J], we release the assumption of equal nucleon masses. This means 
that particle one and two can have different masses which leads to additional independent 
matrix elements. Therefore, one has to consider eight instead of six independent partial 
wave or helicity matrix elements for scattering of positive energy states. For completeness we 
present in the following the complete sets of matrix elements for the various OBE amplitudes. 



We follow the notation of ref . [43j . The potential expressions are presented in terms of helicity 
states. First, the expressions for the scalar mesons, the a meson and the 5 meson, are given 
here 

<+ + \V s J \ + +> = C s <+ + \<j) 3 \ + + > (if + if) (Al) 

< + + \V s J \--> = C s <+ + -->(/«- if) (A2) 

< + " \V S J \ + - > = C s < + - + - > (if + /j 0) ) (A3) 

< + -\V/\- + > = C s < + -\<p s \-+>(lf -if) (A4) 
<+ + \V/\ + -> = -C s < + + \(f> s \ + -> if (A5) 
<+ + \V s J \-+> = -C s < + + \4> s \ + -> if (A6) 

< + - \V s j \ + +> = -C s < + - \4> s \ + +> if (A7) 

< - + IK J | + + > = -C s < + - \4>s\ + +> if, (A8) 

where one has 

Cs = 7rg 2 NNs (A9) 

and 

< a;a< W a,a 2 > = (i - - (aio) 

e x ei e 2 e 2 

with 6i = E* + m*. The integrals over the Legendre polynomials if — if are those given 



in appendix B of 



43|]. 



Secondly, the expression for pseudovector mesons, the n meson and the rj meson, are 
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written as 
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(All) 
(A12) 
(A13) 
(A14) 
(A15) 
(A16) 
(A17) 
(A18) 



(A19) 



(A20) 



with scaling mass M. In eq. (1A20I) the Blankenbecler sugar or Thomas approximation is 
used, i.e. the exchanged energy transfer between the two nucleons is restricted to zero. 
Therefore, the four-momentum transfer is (p' — p)^ = (0, p' — p). This approximation is 
later on also applied for the vector mesons. 

Finally, the vector mesons, the uj meson and the p meson, are treated. The vector-meson 
exchange potential V v consists of three terms: the vector- vector contribution V vv , the tensor- 
tensor contribution Vu, and the mixed vector-tensor contribution V v t- The vector- vector part 
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can be written as 



< + + \v v J v \ + + > = c4< + + |0 o | + + >(/i 1) + /i o) ) 

+ < + + |0v| + + > (I? - 3J< 0) )] (A21) 

< + + \K J V \ --> = C vv [< + + |0 O | - - > - /f) 

+ < + + |0 V | - - > (J< 1} + 3/f) (A22) 

< + - iKtl + - > = C vv [< + - |0 O | + -> + < + - |0v| + - >](/f + lf ] ) (A23) 

< + - IKil - + > = C vv [< + - |0 O | - + > + < + - |0 V | - + >](/f - /i 0) ) (A24) 

< + + iKil + - > = -CU< + + l^o I + -> + <+ + |0 V | + - >]/f (A25) 

< + + \V V J V \ - + > = -C7 W [< + + |0 O | + -> + <+ + |0 V | + - >\lf (A26) 

< + - \V£\ + + > = -C vv [< + - |0 O | + + > + < + - |0v| + + >]/f (A27) 

< - + iKil + + > = -C w [< + - |0o| + + > + < + - |0v | + + >}I?\ (A28) 

where one has 

C vv = ng 2 NNv , (A29) 
< A 1 A 2 |0o|AiA 2 > = (H ; )(H ; ), (A30) 

and 

^ I'i'U ix x ^ 2\' lP ' 2\ lP 2\> 2 p> 2\ 2 p 

< A 1 A 2 |0 v |AiA 2 > = -( — -, 1 )( — 1 ). (A31) 

€ 1 Ci e 2 62 



28 



The tensor-tensor part is 

< + + \vd\ + + > = c tt [< + + 101*1 + + > + /f) 

+ < + + |0i,| + +>(/J 4) + /j 1) ) (A32) 

+ < + + |0 CT *| + +>(/j 1) -3/f)] 

<+ + \Vj\--> = C tt [< + + |0 lt | - - > (/« - if) 

+ <+ + -->(/i 4) -/J 1} ) (A33) 

+ < + + |*«| - - > (I? + 3/f)] 

< + - \Vd\ + ~ > = C tt [< + - |0i*| + -> + < + - |<U + - >](J< 2) + if) 

+ < + " I0i.l + ->(/i 5) +/i 1} ) (A34) 

< + - \V/ t \ - + > = C tt [< + - |0 lt | - + > + < + - \4> at \ - + >](/f - j<°>) 

+ < + - |0i,| - + >(/f-/J 1) ) (A35) 

< + + \Vd\ + - > = -C tt [[< + + + -> + <+ + \4>at\ + - >}lf 

+ <+ + \<f> w \ + -> if] (A36) 

< + + \V/ t \ - + > = -C tt [[< + + |0u| + -> + <+ + |0 CT *| + - >}lf 

+ <+ + |0i e | + ->/f] (A37) 

< + -\V t i\ + +> = -C tt [[< + - |0 lt | + +> + < + - \<f> at \ + + >}lf 

+ < + - |0i,| + +>/f] (A38) 

< - + \V t i\ + + > = -C tt [[< + - |0i*| + +> + < + - |0 CT *| + + >}lf 

+ < + - |0i,| + +>/f], (A39) 



where one uses 



Cu = x*™ (A40) 



Furthermore, one has 
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with A tt = (ml+m'f)(mZ+m$), B tt = (ml+m[*) 2 +(m*+m'*) 2 +(E[*+E*)(E 2 *+E*)+p 2 +p' 2 , 
D tt = -(E' 2 * + E* + E* + E'*)(ml + m'*), and E tt = -(E' 2 * + E* 2 + E{ + E'*)(m* 2 + m'*). 
Furthermore, the mixed vector-tensor part is given by 
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< + - 




+ + > 


< - + 




+ + > 



c vt \< 



+ -> + < + - 10,1 + - >](/f + if) (A46) 



Kll - + > = C vt [< + - - + > + < + - |0 CT | - + >](jf - J^) (A47) 



r(0)N 



,(3) 



CJ< + 



k\ + ~> + < + + |0 CT | + " >]/j 



(3) 



\ ; ! + -f > = -c vt [< + - |0x| + + > + < + - |0 CT | + + >]/f 
K*l + + > = -Cvt[< + - |0i| + + > + < + - |0,| + + >]/i 3) , 



(A48) 
(A49) 
(A50) 
(A51) 
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where one has 



and 



Cvt = 71 2M ' (A52) 



< A 1 A 2 |0i|AiA 2 >= [A vt + D vt — J, (A53) 

e l e 2 e l £ 2 



< A / 1 A / 2 |0 <J |AiA 2 > = (ml + mi* + m 2 + m' 2 *) < AiA' 2 |0 v |AiA 2 > 

+ [W _ m ^P _ 2A^ )( 2A^ + 2A^ } (A54) 
ei e x e 2 e 2 

2A lP 2A' lP ' 2A 2 p 2A' 2 p\ 

+ \P2 -h 2 )(—— + — — )(— 3~ J 

e i e i e 2 e 2 
with A vt = 2(m'? + m 2 * + + - - £ 2 * - E{ - E%) and D vt = 2(e' 1 + e' 2 + d + e 2 ). For 
the uj meson only the vector- vector part contributes, because ()NNt = for the uj meson. 



APPENDIX B: PARTIAL WAVE DECOMPOSITION 

For a general two-body reaction with four distinct spin-1/2 particles and ignoring anti- 
particles, the number of independent amplitudes is sixteen. Due to parity conservation this 
number is reduced to eight independent amplitudes. We denote a helicity amplitude by 
< A' 1 A 2 |0 J (p',p)|AiA 2 >, where \ and A- are the initial and final helicities, respectively. 
Therefore, we have the following set of amplitudes 





=< + + W\p\p) 


+ + > 


<t>i{p',p) 


=< + + \<P J {p',p)\ 


— > 


<t>l{p',p) 


=< + -\(f) J (p',p)\ 


+ - > 


<t>l{p',p) 


=< + -\<P j (p',p)\ 


- + > 


4>1{p',p) 


=< + + \<f) J (p',p)\ 


+ - > 


<f>i(p',p) 


=< + + \<P j (p',p)\ 


- + > 


<Pi(p\p) 


=< + -\<P j (p',p) 


+ + > 


Pt(p',p) 


=<- + \<P J (p',p)\ + + > 
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To partially decouple this system, it is useful to introduce the following linear combinations 
of helicity amplitudes 



rpj 

1 12,0 ~~ 


ti- 


2 J 


rpj 

J 34,o ~~ 


03 J " 


04 


rpj 

1 56,o — 


05 J - 


06 J 


rpj 

J 78,o _ 


07- 


08 


rpj 

1 12,e — 


01 + 


02 J 


rpj 

J 34,e — 


03 J + 


<t>l 


rpj 

J 56,e _ 


05 + 


06 


^78,e = ^ 


^7+</ 





(B2) 



In solving the coupled scattering equation using the linear combinations of the helicity 
amplitudes, two subsets of coupled integral equations, 



T 12 = V 12 4- / v 12 T 12 4- T/ 56 T 78 

o o I o a ' o o 



Tf = Vf + / V a M T? + V?T? (B3) 



rp56 T/56 _i_ / t^12^-t56 1 t^56t->34 

o *o 1*00 '00 



and 



"■78 _ t/78 1 / T /78ml2 , yUrp7% 
o o ' ' o a 



T] 2 = V e 12 + I V} 2 T} 2 + K 56 Tj 8 



Tf = V? + / V^Tf + V™Tf (B4) 



^->56 _|_ J yl2rp56 _|_ ^56^->34 



Tj 8 = Vj 8 + I V™T} 2 + V P U TJ 8 



emerge. Equation (1B3I) is a coupled spin singlet-triplet state, whereas eq. (1B4|) is coupled 
triplet state. For identical particles the coupled spin singlet-triplet state of eq. (IB3I) decouples 
further into a decoupled singlet state 

T c 12 = V a 12 + / K 12 T C 12 (B5) 
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and a decoupled triplet state 



T; 



34 



r o 1 



y 34^34 
* o J o 



(B6) 



due to 56 V J 



_78 



J _56 J,J 



.78 



Tl = 0. 
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